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Abstract 

We explore in detail the physics potential of a measurement of the tt invariant 
mass distribution. First, we assess the accuracy of the best available predictions for 
this observable and find that in the low invariant mass region, the shape is very well 
predicted and could be used to perform a top mass measurement. Second, we study 
the effects of a heavy s-channel resonance on the tt invariant mass distribution in a 
model independent way. We provide the necessary Monte Carlo tools to perform the 
search and outline a simple three-step analysis. 



1 Introduction 



The top quark is unique among the so-far-discovered matter constituents: it is the only 
fermion whose mass is very close to the scale of electroweak symmetry breaking (EWSB), 
nit ^ v/\^, and therefore it has a naturally strong coupling, Aj ^ 1, to the Higgs boson in 
the standard mo del0 It is because of this large coupling, for instance, that the Higgs mass 
can be predicted via precision measurements, and the Higgs can be copiously produced 
at hadron colliders operating at the TeV scale, via its top-loop mediated interactions to 
gluons. Thanks to its large mass, the top quark has also been exploited in many scenarios 
that go beyond the standard model (BSM). The simplest one is SUSY, where top has the 
important role of triggering EWSB and the historical merit of having escorted the MSSM 
to the LHC era by allowing the Higgs mass to survive to the LEP bounds. 

Top has also been the source of inspiration for many alternative mechanisms of EWSB, 
where it has a direct or indirect role in the boson as well as fermion mass generation: 
Topcolor [H El El H] , top see-saw [HIE], [7] and other refinements [H]. More recently, in other 
theoretical constructions that aim at the stabilization of the Higgs mass, the presence of 
the top quark requires a mechanism that compensates for its large (negative) radiative 
corrections to the Higgs mass. In general such mechanisms entail the existence of new 
particles, i.e., top partners, that can be scalars, like in SUSY, as well as fermions, like in 
the Little Higgs [SlIinilllliniiniElIISlDE] or models with extra dimensions [HI [IB]. 
In addition, in many such models gauge interactions exist whose coupling with the third 
generation quarks and in particular to the top quark are enhanced. These include Kaluza- 
Klein (KK) excitations of the graviton fi9[ [20l [2T] as well as the weak [22] and the strong 
gauge bosons [23l[2ll[25l|26l[T9l[27l[2Hl[29l[30| which couple to top quarks. Such particles 
could show up as resonances in the pp ^ X ti production channel and not in other 
channels, like di-jets or di-leptons, due to their small couplings to light particles. 

If on the one hand, the fact that many rich and theoretically motivated models predict 
new physics in connection with the top quark provides a strong motivation for detailed 
experimental investigations, on the other it makes difficult to perform the analyses cor- 
responding to all the suggested scenarios. A way to avoid such an "explosion" of models 
is a bottom-up approach, where first a physical observable, which carries some potential 
for BSM studies, is identified and then the effects due to the existence of new physics are 
systematically explored. The aim of this paper is to present such a model independent ap- 
proach to the discovery and identification of new physics in rriu, the top-antitop invariant 
mass distribution. We do not focus on specific models, instead we assume the existence 
of heavy particles, whose masses and quantum numbers are unknown but are such that 
interactions with the top quark are privileged, a feature common to many scenarios where 
the top quark plays a role in the EWSB, as discussed above. To simulate the physics as- 
sociated with the new states, we have developed a dedicated "model" in the multipurpose 
Monte Carlo generator MadGraph/MadEvent [^[52l[33]FI 

^Another numerical coincidence is also often mentioned, namely mt ~ mw + mz, with somewhat less 
inspirational effects. 

^The Monte Carlo and a wide collection of parton-level data samples (Les Houches format) suitable 
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The strategies and the difficulties associated to the accuracy which will be needed for 
the reconstruction of rriu have been the subject of several investigations [Ml [35] and are 
left to dedicated experimental studies. However, we briefly comment on them in Appendix 
A. 

The paper is organized as follows. In Section [21 we perform an analysis of the theoretical 
uncertainties on the QCD predictions for mu, in the low- as well as high-mass regions. We 
find that the shape of the mu distribution has very little theoretical uncertainties (in 
contrast to its normalization, i.e., the total cross section) and we argue it could provide a 
new handle on the top mass determination. 

In Sections [3l [H and [5] we explore a three-step analysis. In Section [3] the effects on 
rritt induced by new heavy resonances in the s-channel, pp X ti, from the presence 
of simple peaks, to non-trivial patterns arising from interference between signal and SM 
background. In the following section we assume that a resonance is found and study how 
the spin structure of a resonance affects the angular distributions of the top and anti-top 
quarks. In Section [S| we discuss how spin-correlations are affected by nature of the coupling 
of the resonances to SM particles, by simulating the full matrix elements including the 
decay, pp ^ X ^ ti ^ 6f . We leave to Section [6] the discussion of the results and our 
conclusions. 



2 SM theoretical predictions at NLO 

In this section we study the theoretical uncertainties on the available predictions for the 
invariant mass spectrum and their dependence on the top mass. We mainly focus on the 
LHC and refer to the results of Ref. [36] for the Tevatron. We start by considering the 
invariant mass spectrum of the tt pair calculated up to next-to- leading order (NLO), as 
implemented in MCFM [37]. We use the CTEQ6M PDF-set [38] and do not apply cuts 
on the final state particles. Here and in the following we always assume that the invariant 
mass can be fully reconstructed, see also the Appendix. In Fig. [1] results at the LHC are 
plotted for three different top quark masses with uncertainty bands associated to the PDF 
errors and renormalization and factorization scale variations. 

The PDF uncertainty is estimated by running the 41 members of CTEQ6 PDF set, 
with the scales set equal to fi^ = fip = rrit, and found to be about ±3.2%. The scale 
uncertainty is obtained by varying independently the renormalization and factorization 
scales in the region between fiji = fip = fnt/2 and fJ^R = fJ^r = '2mt. The associated total 
scale uncertainty at NLO is about ±13%. Thus, the theoretical errors at the LHC are 
completely dominated by the scale uncertainty. This is contrast to the Tevatron where 
scale and PDF errors are comparable, of the order 6% [36] . 

NLL resummed calculations suggest that the dependence on the scales could go down 
to ±6% [391 [10], however in that analysis the scales were not varied independently. A more 



for further experimental analysis on the scenarios presented in this work, are available on the MadGraph 



servers, e.g., http://madgraph.phys.ucl.ac.be 
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Figure 1: Scale (dashed) and PDF (dotted) uncertainties in the tt invariant mass spectrum 
for top masses (a) rrit = 165 GeV, (b) nit = 170 GeV and (c) nit = 175 GeV at NLO for 
the LHC using the CTEQ6M pdf set. 
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Figure 2: The MC@NLO mu distribution compared with the LO (blue) and NLO (red) 
fixed order predictions. The distributions are normahzed to the MC@NLO cross section. 
We set mt = 170 GeV, include no cuts, and use CTEQ6M for the NLO and MC@NLO, 
and CTEQ6L1 for the LO calculations. 



recent study suggests that changing the scales independently might increase the error to a 
size similar to that estimated through the NLO fixed order calculation |41j . 

Next we compare the NLO shapes for the invariant mass distribution with those ob- 
tained at LO and MC@NLO [12], Fig. H both at the Tevatron and the LHC. We find 
that the differences are minimal at the LHC, and well within the uncertainty bands of 
the theoretical errors on the NLO cross section. On the other hand, at the Tevatron the 
differences between LO and NLO, fixed and dynamic renormalization and factorization 
scales are larger. 

In the high invariant mass region for the LHC, Fig. [3l the LO approximation starts to 
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Figure 3: Scale (dashed) and PDF (dotted) uncertainties in the tt invariant mass spectrum 
for mt = 170 GeV at NLO in QCD, with the CTEQ6M PDF-sets. Also plotted are the 
LO distribution [light dash-dotted), the LO including NLO Electro- Weak corrections [dark 
dash-dotted) with CTEQ6L1 PDF-set. The LO distribution is normalized to the NLO 
total cross section. 



deviate from the NLO order and clearly underestimates the NLO distribution (note that 
curves here are normalized to the total cross section at NLO). Also, as expected, the PDF 
uncertainties start to increase and dominate the theoretical errors as the most important 
contributions come from the large x region. Next-to-leading order electroweak corrections 
to the LO distribution are also included in this figure [l3l HI]. Their effect is to decrease 
the cross section by a few percent for invariant masses below 1000 GeV and up to 15% 
for invariant masses around 4 TeV (the Higgs mass dependence is mild). This means that 
EW effects on this distribution are negligible compared to the current PDF uncertainties 
and give only a minor deviation from the LO curve. 

We conclude this section by mentioning the other sources of potentially large uncertain- 
ties in the determination of the tt invariant mass. The first is related to its reconstruction 
from the decay products. In general the uncertainty on the rriu distribution will depend 
on the final state signature (fully-hadronic, single- lepton and double-lepton final states), 
which determine the reconstruction technique and, more importantly, on the detector ef- 
ficiencies and resolutions. For completeness we briefly discuss the current proposals for 
reconstruction in the various decay channels in the Appendix. The second is due to both 
QCD backgrounds, i.e. multi-jet, W,Z +}ets and iyiy-|-jets, and top backgrounds, i.e. 
single-top and ti itself as coming from a final state different signature than the one con- 
sidered. While the QCD backgrounds at the Tevatron are severe but very well studied, it 
has been shown that at the LHC their impact at low ti invariant mass is negligible when 
at least one lepton is present in the final state [15]. In the high invariant mass tail, some 
QCD backgrounds, and in particular W+ one or two jets, become important due to the 
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Figure 4: The ti production cross section as a function of the top quark mass rrit including 
scale dependence at the Tevatron (a) and the LHC (b). 



fact that the tops are highly boosted and can give rise to single jet-like topologies when 
they decay hadronically. The interested reader can find a detailed study for the single 
lepton final state signature in Refs. [351 HS]- 

2.1 Top quark mass dependence 

As can be clearly seen from Fig. [1], the normalization, as well as the shape of the ti 
invariant mass distribution depends on the mass of the top quark. It is then natural to 
wonder whether such a rather strong dependence could provide another way to determine 
the mass of the top quark. The aim of this subsection is to provide a quantitative answer, 
based only on the theoretical uncertainties. 

In Fig. m we have plotted the ti production cross section a as a function of the top 
quark mass at the Tevatron (a) and the LHC (b). The scale uncertainties, even at the 
NLO, are rather large. Neglecting non-linear terms, a fit to the central curve gives 

Am^/mt ~ 0.2Aa/(T + 0.03 (LHC). (1) 

This equation relates the relative uncertainty on the measurement of the ti cross section 
to the relative uncertainty on the top quark mass: the Aa/a term represents the slope 
and the constant term the horizontal spread, i.e., the scale uncertainty, of the curves in 
Fig. m This means that a measurement of the cross section with an uncertainty of 5% 
would lead to a 0.2 x 5% + 0.03 = 4% uncertainty of the top quark mass, the error being 
mainly associated with scale variations. At the Tevatron the situation is slightly different. 
The scale dependence is milder, 

Amt/mt ~ 0.2Aa/a + 0.015 (Tevatron) , (2) 
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Figure 5: The normalized tt production cross section as a function of the tt invariant mass, 
rritt, for the Tevatron (a) and the LHC (b). Sohd hues from left to right are for a top quark 
mass of nit = 160, . . . , 180 GeV in steps of 5 GeV, respectively. The bands spanned by the 
red lines show the scale uncertainties. 



and known to be reduced at NLL [36], but the PDF errors, which are not included in the 
plot, are not negligible and are found to be of a similar size [36] . 

We can therefore conclude that the accuracy of an independent extraction of the top 
mass from a measurement of the cross section is limited by the NLO theoretical uncer- 
tainties. Recent work suggests that inclusion of NNLO corrections could reduce the scale 
uncertainties sizably [17]. 

It is therefore worth investigating whether information on the top mass can be extracted 
from some other quantity besides the total cross section. In Fig. [5] the ti invariant mass 

, are plotted for five different top quark masses. 



da 



distributions normalized to unity, 

rrit = 160 . . . 180 GeV in steps of 5 GeV. The bands spanned by the red lines show the 
left-over scale uncertainties which are sizably reduced compared to Fig. [TJ We find that the 
shape of the mu distribution is quite insensitive to theoretical uncertainties, while retaining 
a strong dependence on the top quark mass. It is therefore interesting to consider whether 
the invariant mass distribution could provide an independent measurement of the top quark 
mass. 

One way to quantify to which extent the shape is sensitive to the top mass vs. the 
theoretical uncertainties is to perform an analysis based on the first few moments of the 



normalized tt invariant mass distributions 



This approach has the virtue of being 



simple and systematic. Needless to say, alternative quantities, such as the peak position, 
or more sophisticated techniques, such as Kolmogorov tests, could also be employed. 

In Fig.[6]we present the mean value {mu), standard deviation s, skewness 71 and kurtosis 
72 of the ti invariant mass distributions as a function of the top quark mass rrit for various 



6 



560 



540 



520 



500 



480 



460 



- 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

<ni^f> [GeV] 




' ' 1 - 


- 

— 
- 


NLO, CTEQ6M, 


LHC 


- 

— 
- 




Amj/iiij ~ 1.2 A<mjj>/<mj.j> 








+ 0.003 ^ffl35^^ 








^^^^ Red: 


l^n = 






^^^^'^ Black: 


f^s = 






Blue: 


A'r = 


2mj — 




Dashed 


= 


m/s : 




Solid: 


Mr = 






Dotted: 


f^F = 


2m, 








. . \ '- 



160 



165 



170 175 180 

top quark mass, [GeV] 



(a) 



1.35 



1.30 



1.25 



1.20 



I ' ' ' ' I ' ' ' ' I 
Skewness 

NLO, CTEQ6M, LHC 



Am/m, ~ 0.7 Ay^ + 0.005 




Red: 


A'r = 


m,/2 


Black: 


f^s = 




Blue: 


= 


2m, 


Dashed: 


Mr = 


m,/2 


Solid: 


Mp = 


m. 


Dotted: 


Mp = 


2m, 



160 



165 



170 175 180 

top quark mass, m, [GeV] 



150 



145 



140 



135 



130 



125 



120 



Standard Deviation [GeV] 
NLO, CTEQ6M, LHC 



m,/2 



Red: = 
Black: = m. 
Blue: /ip - 2m, 



Dashed: fi^ = m,/2 
Solid: fJ-p = m, 
Dotted: /ip — 2m, 



160 



165 



170 175 180 

top quark mass, m, [GeV] 



(b) 



I ' ' ' ' I 

Kurtosis 7g 

NLO, CTEQ6M, LHC 
Am,/m, ~ 0.2 Ay^ + 0.005 



1.4 ■ 



1.0 — 




Red: 


Mr = 


m,/2 


Black: 


Mr = 


m. 


Blue: 


Mr = 


2m, 


Dashed: 


Mr = 


mt/2 


Solid: 


Mr = 


m. 


Dotted: 


Mr = 


2m, 



160 



165 



170 175 180 

top quark mass, m, [GeV] 



(c) (d) 

Figure 6: The average value (a), standard deviation (b), skewness (c) and kurtosis (d) of 
the ti invariant mass distribution as a function of the top quark mass rrit including the 
scale dependence at the LHC. 



scales at the LHC. We remind that the skewness (kurtosis) is a dimensionless quantity that 
gives a measure of asymmetry (peakedness) of a distribution. These quantities are defined 

as 
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respectively. The central moments fin are defined as 
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In our analysis we focus on the low invariant mass region and therefore we have limited 
the rriu integrals to rriu < TricntoS = 1 TeV. The aim of this cut is just to mimick an 
experimental analysis where the precision on the higher moments would be limited by the 
statistics. Since our purpose is only for illustration we do not consider these effects further. 
However, we stress that our numerical results do retain a significant dependence on this 
cutoff. 

Due to the small scale uncertainty and the strong linear correlation, the mean of 
the tt invariant mass distribution, (m^j), Fig. Et^a), appears to be an excellent estima- 
tor of the top quark mass. From the experimental point of view, one can also hope 
for smaller uncertainties than those associated to the measurement of total cross section. 
In fact, to measure the mean, many systematics, such as those coming from luminos- 
ity or tagging efficiencies, are much less important. A fit to the mean value shows that 
Amt/rrit ~ 1.2A{'mtt) / {rriu) + 0.003. So, for instance, if the mean value is measured with 
a 1% uncertainty, the uncertainty of the top quark mass is only 1.3%, including the scale 
uncertainties. 

The standard deviation. Fig. [6](b), is almost constant and therefore is not suitable for 
a top quark mass measurement. In Figs. [6](c) and[6](d) the skewness and the kurtosis for 
the tt invariant mass are plotted, respectively. Also here, the scale uncertainty is reduced, 
while still slightly larger than for the mean value. Fig. [6]^a). The slopes of the lines are 
promising, in particular for the kurtosis, which means that a relatively large experimental 
error on the measurement of the kurtosis leads to an only small error on the top mass 
measurement. 

At this point, we have to stress that the above simple analysis does not include neither 
statistical nor systematics effects in the data, which should also be carefully considered. 
In particular, the higher monents such as the skewness and the kurtosis are more sensitive 
to the tail of the rriu distribution then the lower moments and therefore more sensitive to 
statistical and systematic effects that affect more strongly this tail. Eventually, the final 
uncertainty on the top quark mass will depend on how well the above quantities can be 
measured. It is plausible to expect that a combined analysis based on the above quantities 
might lead to an even smaller uncertainty for the top quark mass. 

For completeness we show the same analysis performed at the Tevatron energies, see 
Fig. [71 Also in this case, we have used a fixed order NLO calculation to estimate the scale 
uncertainties. However, as we have already mentioned, at the Tevatron the tt pairs are 
produced almost at threshold, hence a resummed calculation which predicts a smaller scale 
uncertainty, is preferred [M]- 

The Tevatron results (Fig. [7]) are similar to those obtained in the LHC study, but the 
reduction in the scale uncertainties by analyzing the (higher) moments is smaller compared 
to the LHC. The first moment, i.e., the mean value, is probably the best estimator for 
the top quark mass among all the moments, due to its small constant value of 0.004, 
and the reasonably good proportionality factor of 1.2. The higher moments are more 
sensitive to statistical fluctuations and might be less suitable with a limited sample. The 
lack of events in the higher invariant mass regions might give rise to larger errors for the 
skewness and worse the kurtosis. In these plots we restrict the ti invariant mass to below 
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Figure 7: The average value (a), standard deviation (b), skewness (c) and kurtosis (d) of 
the ti invariant mass distribution as a function of the top quark mass rrit including the scale 
dependence at the Tevatron. For the skewness and kurtosis we restricted the integration 
region in Eq. H] to m^f < 600 GeV. 



^cutoff = 600 GeV. We mention that even though using the fixed order NLO calculation 
we have overestimated the scale uncertainties, we have neglected the PDF errors which at 
the Tevatron can reach the 6-7% [36] and errors coming from the reconstruction of the tt 
invariant mass from the (anti-)top quark decay products, see the Appendix. 

Finally, we comment about the definition of the top quark mass. As in the more stan- 
dard top mass measurements where the top mass is reconstructed from its decay products, 
also the top quark invariant mass is sensitive to extra radiation and to non-perturbative 
effects due to confinement, typical of a pole mass [H]. In this respect the same issues and 
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problems in associating a theoretically well-defined mass to the measurement remain. This 
is at variance with an extraction of the top mass from a cross section measurement which 
can be directly related to a short distance mass and does not suffer from the same non- 
perturbative or extra radiation effects. We stress, however, that mu is at least twice the 
top quark mass, which would decrease the relative impact of the ambiguities due to extra 
radiation. In addition, typical combinatorial systematics associated to the assignment of 
the jets to the "right" tops, are absent for m^f. In this respect, more experimental work 
on the systematics affecting such a measurement would be certainly welcome. 

3 Effects from BSM resonances 

In this section we investigate the effects of (model-independent) new resonances on the 
ti invariant mass spectrum. All the numerical results presented here have been obtained 
with MadGraph/MadEvent, through the implementation of a dedicated "model", topBSM, 
which is publicly accessible on the MadGraph servers for on-line event generation and for 
download jfl 

topBSM offers the possibility of studying a wide range of new physics resonances and 
efficiently exploits the flexibility and the possibilities of MadGraph: 

• SM effects are consistently included, i.e., possibly non-trivial interference effects be- 
tween new resonances and the ti background are taken into account. As it will be 
shown in the following, in some cases such effects can be important and might lead 
to very distinctive signatures (cf. the case of the peak-dip structure arising in rriu 
due to the presence of a (pseudo-) scalar state). In general, they should be always 
included. 

• The full matrix elements 2 — 6 including the decays of the top quarks can be 
generated, which is crucial for spin correlation studies. 

• The generated events can be automatically interfaced to parton showers programs, 
such as Pythia [19] or Herwig [50], to shower and hadronize the events after which 
these events can be processed by a detector simulation for full experimental analyses. 

We have considered s-channel spin-0, spin-1 and spin-2 resonances, of different color 
and CP parity, as listed in Tabled! The parameters related to each resonance are simply 
the mass, the width and the relevant values of the couplings to standard model particles 
which enter in the production process (to the partons and to the top quark). 

3.1 Spin-0 resonances 

The first resonances we discuss are spin-0 particles. We distinguish between color singlet 
(0) and color octet (5*°), as well as parity even (scalar) and odd (pseudo-scalar) spin-0 

•^Technical documentation on how to use the model can be found at 
http : / / cp3wks05 . f ynu . ucl . ac . be/twiki/bin/ view/ Sof tware/TopBSM[ 
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Table 1: The BSM particles included in the topBSM "model". 

particles. 

3.1.1 Color singlet 

Let us start by considering a color singlet (pseudo-) scalar boson contributing to the tt 
process gg (0 -^^tt. The Feynman diagram for this loop induced process is depicted in 
Fig. [HI The spin-0 coupling strength to quarks, 

g^q = '^1^ ^ ^2— 75, (5) 

V V 

is proportional to the quark mass m^. In analogy with the SM, v is the spin-0 field vacuum 
expectation value and ai and 02 are real proportionality factors for the parity even and 
odd spin-0 particles, respectively. For the SM Higgs boson ai = 1 and 02 = 0, while for a 
pure pseudo-scalar oi = and 02 is non-zero. 

We do not include scalar production by (anti-)quark annihilation, qq 0, because for 
this cross section to be sizeable compared to the loop induced gluon fusion process, the 
branching ratio for the scalar to ti has to be small and can be neglected. 

Since we are interested in scalars with strong couplings to the top quark, we neglect 
all particles in the loop of Fig. [8] except for the most heavy quark, i.e., the top quark. 
If the mass of the spin-0 boson is larger than twice the mass of the top quark, the loop- 
induced gluon-gluon-(pseudo-)scalar coupling develops an imaginary part, which leads to a 
peak-dip structure for the interference terms between the QCD background and the signal 

The possibility to detect a signal in the ti invariant mass depends on the width of the 
spin-0 resonance. In general, a scalar particle couples also to the electroweak bosons. In 
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Figure 8: Feynman diagram for the (pseudo-) scalar contribution to tt production. 



the SM the decay rate to VF, Z is much larger than the decay rate to tt, and therefore 
the tt channel is suppressed. Moreover, the presence of a destructive interference between 
the signal and the QCD background and the relatively large width of the scalar makes 
detection very difficult. An enhanced coupling to top would not help much because the 
improvement in the branching ratio would be compensated by an increase of the total 
width. In conclusion, there is little hope to see a SM-like scalar by looking at the tt 
invariant mass spectrum, even if the coupling to top quark were (much) larger than in the 
SM. 

On the other hand, the case of a pseudo-scalar or a 'boson-phobic' scalar resonance 
that does not couple to the heavy vector bosons is more promising. For such a state, the 
branching ratio to tt can be taken unity, BR(0 ^ tt) = 1, i.e., the total width of the scalar 
spin-0 resonance is equal to the SM partial width to ti. SUSY models with this feature can 
be constructed [59]. The smaller widths of the pseudo-scalar and the boson-phobic scalar 
give a narrow resonance peak in the ti invariant mass spectrum. The interference between 
the signal and the QCD tt production leads to a dip in tt production at an invariant 
mass just above the mass of the spin-0 particle. In this case the signal together with the 
interference terms sum to the characteristic peak-dip structure. Fig. [91 

The dot-dashed line in Fig. [9] shows the effect of a 400 GeV color singlet spin-0 particle 
on the tt invariant mass spectrum with couplings ai = 1, a2 = and ai = 0, a2 = 1 for the 
left and the right plots, respectively. Comparing with the QCD ti production, the dark 
solid line, a peak-dip structure is visible when the spin-0 particle is a pseudo-scalar, ai = 
and a2 = 1. In the case where it is a scalar, ai = 1 and 02 = 0, there is only a peak and a 
very small dip. 

If the coupling to the top quark is enhanced, the peak as well as the dip becomes 
broader due to the larger decay width. The peak increases in the case the spin-0 is a 
scalar, but remains the same for the pseudo-scalar. The dashed line shows the effect of 
enhancing the ttH coupling by a factor of two. If the coupling to the top quarks is taken 
even larger, the increasing width of the (pseudo-) scalar starts to dominate the effects on 
the invariant mass. This results in the disappearance of the dip, as shown by the light solid 
line in Fig. [9l 

In the case where the coupling to the top is smaller than in the SM, the peak of the 
scalar gets smaller and the dip completely disappears. The effect of varying the coupling 
for the pseudo-scalar are much smaller. Even if the coupling to top quarks is reduced by a 
factor of two, ai = and 02 = 0.5, a very clear peak-dip structure is still visible, as shown 
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Figure 9: Invariant ti mass spectrum for the boson-phobic scalar {left) and pseudo-scalar 
{right). Bottom: The interesting region with finer binning. Different colors represent 
different coupling strength of the Higgs to top quarks: dot-dashed for the standard model 
coupling and dotted, dashed and light solid for 0.5, 2 and 4 times the standard model 
coupling strength, respectively. Dark solid is QCD tt production, i.e., without the Higgs 
signal. All plots were produced using the CTEQ6L1 pdf set with fip/^ = fip = 400 GeV. No 
acceptance cuts are applied. 



by the dotted line in Fig. [9l 
3.1.2 Color octet 

The case of a color octet resonance is very similar. Here we shall study scalar and a 
pseudo-scalar Sj color octets, similar to those introduced in Refs. [Ml ES] • In these models 
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Figure 10: Invariant ti spectrum for pp ti including a s-channel 5"^ scalar color octet 
(a) and a pseudo-scalar scalar Sj color octet (b) with masses 777,50 = nigo = 400 GeV. 
Dark solid line is QCD tt production, dash-dotted line is with standard coupling between 
the scalar and ti, dashed, light solid and dotted the coupling is 2, 4 and 0.5 times as large, 
respectively. All plots were produced using the CTEQ6L1 pdf set with fiR = (j^f = 400 
GeV. No acceptance cuts are applied. 



the (pseudo-) scalar color octet couples only to quarks, with the same SM coupling but for 
the color ^ 

9so,qq = Vui-^^, and ^^0^^ = rju^i^-^T^j, (6) 

where rju is a coupling proportionality factor and of order 1. The production and decay 
mechanism for the (pseudo-) scalar color octet are similar to the 'peak-dip' color singlets, 
i.e., the resonance is produced through a top quark loop by gluon-gluon fusion, and the 
decay is mainly to top quarks. We find that compared to the 'peak-dip' color singlet 
the 'signal' cross section is 5/72 times smaller, i.e., o-{gg ^ S'}^^ ^ ti) = ^o'{gg 
H ^ ti), the interference between signal and background is 5/12 times smaller and the 
width of the (pseudo-) scalar color octet is 6 times smaller than the width of the 'peak-dip' 
(pseudo-) scalar color singlet. In Fig. [TOlfa) the ti invariant mass is plotted in a model with 
a color octet scalar of a mass of 400 GeV, and in Fig. [TOlfb) for the pseudo-scalar. As 
expected, the results are very similar to the color singlet, see the lower plots of Fig. [91 
The same 'peak-dip' structure is also present for the color octet but it is more pronounced. 
This is mainly due to the smaller width of the (pseudo-) scalar color octet. 

We conclude this section by mentioning that pseudo-scalar singlet and octet resonances 
could also arise from bound states of meta-stable gluinos in split SUSY scenarios [60]. Also 
in this case gg ^ (p would be the dominant production channel and it could be described 
within the same framework. 
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Figure 11: Invariant tt spectrum for "pj) tt including a s-channel Z' color singlet vector 
boson and color octet (axial) vector bosons with masses mx = 2000 GeV that couples 
with standard model strength to quarks. Solid QCD tt production, dotdashed with a color 
singlet [Z'), dotted with a color octet axial vector (axigluon g\), dashed with a color octet 
vector boson (KK gluon/coloron Qy). All plots were produced using the CTEQ6L1 pdf set 
with = Hp = 2000 GeV. No cuts were applied in making any of the plots. 

3.2 Spin-1 resonances 

In this section we discuss a spin-1 resonance produced by qq annihilation. This resonance 
can either be a color singlet or a color octet. For the color octet case we distinguish between 
a vector and an axial-vector. Although both the vector and the axial-vector interfere with 
the QCD ti production, only the vector shows interference effects in the ti invariant mass 
spectrum. 

Including an s-channel color singlet vector boson (a "model-independent" Z') in the ti 
production process gives a simple peak in the invariant mass spectrum as can be seen from 
the dot-dashed line in Fig. [TTl The precise width and height of the peak depends on the 
model parameters in the model for the Z' . As a benchmark we show a Z' vector boson 
with mass m^/ = 2 TeV that couples with the same strength to fermions as a standard 
model Z boson. The interference effects with the SM Z boson can be neglected in the ti 
channel, so the peak is independent of the parity of the coupling. 

In general, for the color octet spin-1 particles the interference with the SM ti production 
cannot be neglected. Two cases are to be considered: a color octet vector particle {e.g., a 
KK gluon [nS] or coloron and an axial-vector particle {e.g., an axigluon [HH E21 EZ] ) • 
It is natural to assume a coupling strength equal to the strong (QCD) coupling for their 
coupling to quarks. 

In Fig. [11] the effects of a color octet spin-1 particle on the ti invariant mass spectrum 
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are presented. The interference effects of the axial vector {dotted line) with the QCD tt 
production does not change the shape of the tt invariant mass spectrum. Hence the effects 
of the color octet axial vector and the color singlet are very similar, apart from the size 
due to the different coupling constant. 

The interference of the color octet vector particle with the QCD tt production does 
effect the tt invariant mass distribution. There is negative interference in the invariant 
mass region below the resonance mass and positive interference for heavier invariant masses. 
This slightly changes the shape of the peak as can be seen from Fig. [TTl Other quantities, 
such as the charge asymmetry between the top and the anti-top quarks could be more 
sensitive to axial vectors [63] and could help their discovery at higher invariant masses. 

3.3 Spin-2 resonances 

The interactions between spin-2 particles, or gravitons, and ordinary matter is in general 
Planck suppressed, which makes it impossible to see effects of the gravitons at TeV ener- 
gies. There are, however, models with extra dimensions where the contributions from the 
gravitons might be large enough to make a discovery at the LHC In this case a model- 
independent approach is not really appropriate. Instead we consider two scenarios that 
have distinct signals in the ti invariant mass. First the ADD model [171 El] , where the 
effect of a the large number of graviton KK states contributing to a cross sections could be 
important and, secondly, the RS model [T8l [65] where only a limited number of KK modes 
contribute, but the coupling constant itself is enhanced by a large "warp" factor. 

In the so-called ADD models [TJl [M] all SM fields are confined to a four-dimensional 
brane, letting only gravitons propagate through the bulk. The extra n bulk dimensions are 
compactified on a n-torus with a radius R. If the radius R is large enough, (of the order 
of 0.1 mm for 2 extra dimensions) the (4 + n) dimensional Planck scale can be as small as 
the TeV scale. 

Due to the fact that the radius of the extra dimensions is large, the graviton KK 
states can be almost degenerate in mass. So, although all graviton couplings are Planck 
suppressed, the sum of all the KK states can contribute significantly to the ti invariant 
mass spectrum. All states are summed up to the cutoff scale Ms, defined by A^i?" = 
8-7r(47r)"/^r(?i/2)M^"^^, where A is related to the four dimensional Newton's constant 
A = v^TBvrG/v- 

The effect of this tower of graviton states on the tt invariant mass distribution is plotted 
in Fig. [12] in the case of 3 extra dimensions and for 4 different cutoff scales. Due to the sum 
over all nearly degenerate resonances, there is no single resonance peak in the invariant 
mass distribution. It is also clear that the distribution is only valid well below the cutoff 
scale Ms, otherwise unitarity violating effects become sizable. 

In the so-called RS model [181 ES] there is one extra dimension postulated that is 
compactified to a S^/Z2 orbifold. There are two branes on specific points of the orbifold: 
a "Planck" brane at = and a "TeV" brane at = vr where the physical SM fields are 
confined. The bulk space is warped in such a way that the (reduced) Planck mass is warped 
down on the "TeV" brane to A = Mpie~^'^^. The gauge hierarchy problem (A = 0{1 TeV)) 



16 



100 



do-(pp -> (G ^) tt)/dm, 
LO, CTEQ6L1, LHC 



■u 



[pb/20 GeV] 



80 



Ms=800 GeV 

Mg=900 GeV 

Mg=1100 GeV 



M_=1300 GeV 



60 — 




SM 



40 — 



20 







400 



600 



800 1000 



tt invariant mass [GeV] 



Figure 12: Invariant ti spectrum for pp — > ti including s-channel gravitons. The distribu- 
tions show the effect of the almost degenerate tower of KK gravitons in the ADD model 
with n = 3 extra dimensions and, from top to bottom, with a cutoff scale Ms = 800, 900, 
1100 and 1300 GeV. The bottom line are contributions from SM only. We used CTEQ6L1 
and set the scales to fj,R = fip = rrit. 

is now solved with only a minor fine-tuning of nR ^ 12. After KK compactification of the 
massless graviton field, the coupling constant of KK gravitons with matter is given by the 
inverse of A. 

A prediction in the RS model is that the masses of the KK modes m„ are given by 
'^'^^ where the positive zero's of the Bessel function Ji{x). If one of the 

masses is given, all the others are fixed, which could give rise to a series of resonances in 
the tt invariant mass spectrum. 

In Fig. [12] the effect of a series of KK graviton modes on the tt invariant mass spectrum 
is shown with mi = 600 GeV and for various ratios n/Mpi. The resonances are clearly 
visible over the QCD background. Higher KK states are characterized by larger widths. 

4 Spin information from (anti-)top quark directions 

A useful, yet simple, quantity sensitive to the spin of the intermediate heavy state into a 
tt pair, is the CoUins-Soper angle 6 [66]. This angle is similar to the angle between the top 
quark and the beam direction, but minimizes the dependence on initial state radiation. 
6 is defined as follows. Let pa and ps be the momenta of the incoming hadrons in the 
rest frame of the top-antitop pair. If the transverse momentum of the top-antitop pair is 
non-zero, then p^ and pb are not coUinear. The angle 6 is defined to be the angle between 
the axis that bisects the angle between pa and pb and the top quark momentum in the tt 
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Figure 13: Invariant ti spectrum for pp — > ti including s-channel gravitons. The distribu- 
tion shows the effect of a couple of KK resonances in the RS extra dimensions model. The 
mass of the first KK mode is mi = 600 GeV and the colored lines represent various choices 
for the ratio n/Mp\. We used CTEQ6L1 and set the scales to fipi = fip = f^t- 



rest frame. 



4.1 Standard Model 

The distribution of 6 in the SM is plotted in Fig. [HI Also plotted in the same figure are 
the distributions with cuts on the tt invariant mass spectrum as backgrounds to narrow 
resonances. 

A simple analytic calculation confirms this behavior. The matrix element squared for 
the initial state qq to the SM ti contribution in terms of the Collins-Soper angle cos 6* is 
proportional to 

\M{qq ^ ti)\^ s{l + cos^ 9) + 4m,^(l - cos^ 9), (7) 

where s is the center of mass energy squared, s = {pg +PqY. For the gg initial state we 
have 



\M{gg^ti)\' 
s{7 + 9cos^i 



36mf cos^ 9 



(sc_ + 4m^ cos^ 6*) 
where c+ = 1 + cos^ 9 and c_ = 1 — cos^ 9. 



s c+c^ + 2sm. ( 3c_ + c 



4m^ (^3c? 
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Figure 14: Normalized distribution for cos 6', where 9 is the Colhns-Soper angle, for SM 
production at the LHC Different lines represent different cuts on ti invariant mass. 

4.2 Spin-0 resonances 

In Fig. [TST a) the normalized cross section as a function of cos^ is plotted for a spin-0 
resonance. The distribution is independent of the mass and parity of the resonance. The 
matrix element squared for the spin-0 resonance H is at leading order proportional to 

\M{gglqq ^H^tt}\^r^ {\a^\^ + \a2\^)pt ■ pt - (|aip - lasHm^, (9) 

where pt and pt are the momenta of the top and anti-top quarks, respectively, and ai and 
a2 are the coupling constants, see Eq. [5l for the scalar and pseudo-scalar, respectively. The 
matrix element squared is clearly independent of the angle cos^, which explains the fiat 
distribution. 

4.3 Spin-1 resonances 

For a generic spin-1 resonance Z' the matrix element squared is proportional to 

\Miqq ^ Z' tt)f ~ 2(|ai|^ \aR\^)Pq- Pt Vq' Vt^ ^\aL^\a-Rfpq- Pt Pg ■ Pt+ 

(IclI^ + \aR\^){aLa*R + aRO^pq ■ pg, (10) 

where ai and are the left and right handed part of the couplings of the Z' resonance 
to quarks, i.e., gz'qq ~ ^^l^-^ + o-r^-^ and where Pq and pq are the momenta of the 
incoming quark and anti-quark, respectively. In terms of cos 6 the matrix element squared 
is proportional to 

\M{qq ^Z' ^ ti)\^ ~ {\aL\^ + \aR\^){s - 4m2)(l + cos' 6) + Am^^laL + a^l^ (11) 
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Figure 15: Normalized distributions for cos ^, where Q is the Colhns-Soper angle, for spin-0 
(a), spin-2 (b), vector (c) and axial- vector (d) resonances of mass Mx- All plots were 
produced using the CTEQ6L1 pdf set with \xr = [if = Mx- 



The normalized cos^ distribution is independent of the mass of the resonance for a axial 
vector, ttR = —ttL (see Fig. [TST d)). while for a pure vector resonance the dependence is 
maximal Fig.fTMc). However, for heavy resonances, Mx ^ 800 GeV the difference between 
the curves for the vector and the axial-vector is less then 8% which makes it challenging 
to get any information about parity of the coupling from this distribution. In Ref. [M] 
a similar polar angle has been studied. That polar angle is also sensitive to the chirality 
of the coupling. However, the CoUins-Soper angle used here has the advantage that it 
minimizes the effects from initial state radiation. 
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4.4 Spin-2 resonances 



In the case of the spin-2 resonance Gf^^, both the qq and gg initial states contribute. The 
matrix element squared for the qq initial state is proportional to 

\Miqq Gf,^ ~ s(l - 3cos^^ + 4cos^^) + 16m2cos^^(l - cos^0), (12) 

and for the gg initial state 

\M{gg ^ Gf,^ ^ti)\'^ [s(l + cos^ ^) + 4m^(l - cos^ (1 - cos^ ^). (13) 

The large differences in the distributions for the spin-2 resonances between light compared 
to heavy spin-2 particles, see Fig. [T5](b). is due to the fact that the relatively light spin-2 
particles are mainly produced by gluon fusion, while the very heavy spin-2 particles by 
quark-antiquark annihilation. 



5 Spin correlations in ( ant i-) top- quark decays 

In the standard model, the semi- weak top-quark decay width is rather large F 1.5 GeV > 
Aqcd and top quarks do not form bound hadronic states. At present, we do not have any 
direct measurement of the top width and the formation of top hadrons is not excluded. 
This could happen for example, if Vtb were much smaller than what is predicted in the 
standard model, as discussed in Ref. [67]. Note, however, that even if this were to happen, 
the information on the spin of the top quark would be anyway fully inherited by its decay 
products [68], as spin-flip would occur at time scales of the order nit/AQ^j-), i.e., much 
later than the lifetime of the top quark. In this respect, spin correlation effects are a very 
robust probe of new physics entering in the production cross section. 

For standard model leptonic top decays, the directions of the leptons are 100% corre- 
lated with the polarization of the top quarks. The spin analyzing power of the direction 
of the b quark {W~^ boson) is not as good, around —0.4 (0.4). In hadronic top decays the 
anti-down (or anti-strange) quarks coming from the boson decay have the same full 
spin analyzing power as the lepton. On the other hand, the up (or charm) quarks have a 
spin analyzing power of only —0.3, i.e., the same as the neutrino in leptonic decays. For the 
decay of anti-top quarks or spin-down top quarks, all spin analyzing powers change sign. 
The angular distributions of the two down-type fermions (leptons in leptonic top decays 
or jets coming from down-type quarks in hadronic W decays) give maximal information 
about the spin of the (anti-)top quarks in ti events [HHl 170] . 

In studies on the spin correlations in tt production, two distributions are usually con- 
sidered [ini EI] . First the distribution 

^-^ — = ifl - A cos 6^+ cos 6'_ + b+cos6+ + 6_ cos6'_V (14) 

a a cos O^d cos 0_ 4V / 

where 6^ {6_) is the angle between the t (t) direction in the tt center of momentum frame 
and the ifd) direction in the t (t) rest frame, where ifd) down-type fermion 
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coming from the decay. For a parity conserving ti production mechanism, such 

as QCD, the parameters 6+ and 6_ vanish. In practice, the way to construct these angles 
is first to construct the t and i four-momenta in the laboratory frame. Then perform 
a rotation-free boost from the laboratory frame to the tt center of momentum frame to 
define the t (t) direction in the tt center of momentum frame. Thirdly, boost the down- 
type fermion momenta, i.e., the lepton in leptonic top decays and the down- type quark 
in hadronic W decays, from the tt center of momentum rotation-free to the t and t rest 
frames. If the t and t rest frames are constructed directly by boosting from the laboratory 
frame a Wigner rotation has to be taken into account [TT] . 

Defining the angles 9+ and 9- as described above, corresponds to studying spin corre- 
lations of the ti pair in the helicity basis. 

It is important to stress that in spin correlation studies it is mandatory to reconstruct 
the top and the anti-top quark momenta. In the case of a double leptonic decay, two 
neutrino's are emitted and the full reconstruction of the event becomes non trivial. Im- 
posing kinematic constraints, such as the known top and W masses, a constrained system 
of equations for the neutrino momenta can be set up. In general multiple solutions arise 
and the best solution can be only obtained on a statistical basis [73l HQ], [TH [21], [TS] [TG] [77] . 
In the appendix some of the issues for the reconstruction are discussed in more detail. 
Alternatively, the single leptonic ti decay could be used by letting one of the jets of the 
hadronically decayed (anti-)top quark play the role of the lepton. Ideally, one would like 
to use the jet coming from the down-type quark, because it has the same (maximum) spin 
analyzing power as the lepton. In experiments one cannot easily distinguish between up- 
and down-type quarks jets on event-by-event basis, and the analyzing power gets averaged 
(1 — 0.3)/2 ^ 0.35. However, improvements can be achieved by exploiting the fact that 
down-type quark jets have in general a smaller transverse momentum than the up-type 
quark jets. Using the least energetic (non-6) jet from the hadronic top decay increases the 
spin-analyzing power from 0.35 to approximately 0.5 [TT] [72]. For illustrational purposes 
in the following we assume that the top quark momenta are correctly reconstructed and 
the spin analyzing power is maximal. 

In Fig. [Tn]this distribution is plotted for QCD ti production. 

The differences among the various ti production mechanisms are manifest. In Fig. [T7] 
the distributions are plotted for resonance masses of 800 GeV for the following states: 

• Scalar boson (a), 

• Pseudo-scalar boson (b), 

• Vector boson (c), 

• Axial- vector boson (d), 

• Vector-left boson (e), 

• Vector-right boson (f). 
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(a) 




Fiffure 16: The distribution --; — ^—^ — s— for SM tt production at the LHC, using pdf set 
CTEQ6L1, without applying cuts (a), and for the regions 390 < niu < 410 GeV (b) and 
790 < mtt < 810 GeV (c). 



• Spin-2 boson (g). 

With the vector-left and vector-right we understand a spin-1 vector boson that couples 
only to left- or right-handed fermions, respectively. We choose very narrow resonances by 
taking the width of resonances to be 1% of the mass, i.e., 8 GeV for a mass of 800 GeV. 
We do not include the SM QCD ti production background in these plots. 

In Table [2] the distributions are fitted to Eq. f|T^ and compared with analytic compu- 
tations. For the sake of simplicity, in the analytic computations the off-diagonal elements 
of the spin correlations matrix in the helicity basis are neglected. This means that the 
interference between different top quark spins are not included. In fact, the interference 
effects are negligible and the fitted values agree very well with the analytic computations. 
For completeness we also included the the numbers for a smaller resonance mass, Mx = 
400 GeV, where the effects from the mass of the top play a larger role. 

The second angle, which is commonly considered when studying spin correlations in tt 
production is </>, i.e., the angle between the directions of the and in the t and i rest 
frames, respectively. The distribution 

a a cos 

for this angle is plotted in Fig. [18] for a resonance mass of 800 GeV. 
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(a) (b) 




Figure 17: The distribution ^ cose +d cose ^'-'^ scalar, (b) pseudo-scalar, (c) vector, (d) 
axial-vector, (e) vector-left, (f) vector-right, (g) spin-2. Mx = 800 GeV at the LHC, using 
the pdf set CTEQ6L1. No cuts were applied. 

The distributions for the angle are the same for production through a scalar and 
a vector boson. The distribution for the pseudo-scalar, on the other hand, is completely 
different from the one for the scalar and the vector boson [53]. Also, the angular distribution 
for SM ti production is different from the other production mechanisms. In the case of a 
spin-1 state, the (f) distribution is independent of the type of coupling to the top: it makes 
no difference whether it is pure vector, an axial-vector, a left-handed or right-handed 
couplings. 
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-0.743 





0.022 
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Table 2: -71 — ^-^ — 5- = 7(1 — A cos 6*+ cos 6*- + b+cos9+ + 6_cos6'_). For a top mass 
of 175 GeV. In the analytic calculation of the parameters, the interference between the 
various top spins is neglected. 

6 Conclusions 

Top physics is entering the precision phase at the Tevatron and will be one of the leading 
priorities at the LHC. The importance of the top quark in the quest for the mechanism of 
EWSB and new physics is due to the convergence of two factors. First the LHC will be 
a top factory with tens of millions of top quarks produced in the first years of running at 
the nominal low luminosity, 2 x 10^^ cm~^ s~^, allowing studies at an unprecedented level 
of accuracy. Second, due to its large mass, top is the optimal probe for new physics at 
TeV scales. Many different models, including the SM, predict the existence of heavy states 
that preferably couple to the top quark, and that could affect its SM couplings via loop 
corrections or be directly produced at the LHC. 

Given the large number of the models proposed and their complexity, a "top-down" 
approach, e.g., model parameter scanning, will not be practical, in particular if comparison 
of many different channels and observables at once will be necessary. As an alternative, 
a simpler and more pragmatic "bottom-up" approach could be employed, whereby one 
identifies specific observables which can be developed as tools for discriminating generic 
features of new physics resonances, thus keeping the analysis as model independent as 
possible. 

In this paper we have presented an example on how such a study could be performed 
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Figure 18: Distribution for the angle 0, defined in the text, for different ti production 
mechanisms at the LHC. Dark Solid line is the SM tt production, dotted line is tt production 
through a scalar, dot- dashed line is tt production through a pseudo-scalar, dashed line is 
tt production through a vector (this is independent of the coupling). The light solid line 
is tt production through a graviton. The plots are normalized. The pdf set CTEQ611 is 
used with Mx = 800 GeV and = = 800 GeV. 

for the invariant mass distribution of the ti pair. 

As a first step we have assessed the accuracy of the best theoretical predictions available 
for tt production at hadron colliders. We have found that the shape of the distribution is 
under good theoretical control, especially at low invariant mass values, suggesting also the 
possibility of a precise top mass extraction. 

We have then identified the features of new physics scenarios, namely the existence 
of heavy bosonic resonances of various spin, color and parity, that could show up in the 
m^j distribution, and implemented them in the MadGraph/MadEvent package. The full 
matrix elements, — > X — > tt — > 6/, X being a spin-0, spin-1 or spin-2 particle, particles 
with arbitrary masses, width, color and couplings, have been automatically generated by 
MadGraph. The effects due to the interference with the pp ti SM process are included 
when relevant. 

The strategy to gain information on new physics is then straightforward and consists 
of three successive steps: 

I. The discovery of the resonance (and the determination of its mass and width) which 
could appear as a sharp or broad peak or as a more distinctive peak-dip structure 
in very specific cases. In this measurement the key aspect will be the experimental 
resolution in the rriu reconstruction. 

II. The identification of the spin of the resonance, which can be inferred from the angular 
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distribution of the top and the anti-top. 



III. Information on the couphngs of the resonance to the the top anti-top pair, which can 
be obtained by measuring the spin correlations of the top anti-top pair (for this last 
step the full matrix matrix element 2 ^ 6 is required). 

In conclusion, we have outlined a simple strategy and provided the necessary Monte 
Carlo tools to search for new resonances in tt events. We look forward to more detailed 
experimental analyses. 
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Appendix: reconstruction issues in tt events 

In this appendix we address some of the issues arising in the reconstruction of the ti events 
and in particular their different impact in the three-step analyses proposed above. 

For a generic 2 — > 6 process, where the final state particles masses are known, 16 
independent variables are needed to determine the kinematics of the event: the six final 
state particle three-momenta, {Pi}, the two energy fractions carried by the initial state 
partons, Xa,b, minus the overall momentum conservation which reduces the number of 
independent variables by 4, leads to 6 x 3 + 2 — 4 = 16. 

From the measured angles and energies of the final state particles, together with con- 
straints from W boson and top quark masses, a system of equations can be set up to solve 
for the 16 unknown variables on an event-by-event basis. 

The three decay modes of the tt pairs face each their own challenges for detection 
and reconstruction. Around 44% of the tt pairs decay hadronically, 30% of the tt pairs 
decay single-leptonically and 5% double leptonic (not including tau's) |4ufi These three 
channels offer very different challenges related to the detection of final state particles and 
the reconstruction the (anti-)top quark momenta, which we will now address channel-by- 
channel. 

The fully-hadronic decays have the advantage that in principle the momenta of all the 
final state particles can be determined, leading to 6 x 3 = 18 measurements. Together 
with the four mass constraints the system of equations for the 16 independent variables 

^The remaining 21% are events including decays to tau's, which are not be considered here. 
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is over-constrained. Such constraints can be used in two ways. First they can be used 
to extract information, typically the jet energies, that have bad detector resolution. For 
example, measuring only the angles [9, (p) of the six jets and including all the constraints 
from the top quark and W boson masses would already provide the required 16 independent 
quantities (althought with combinatorics). Alternatevely, the constraints from the masses 
can also be used to solve the combinatorics in reconstructing the W boson and (anti-)top 
quark momenta. Combinatorics, affect each of the three steps in the analysis proposed in 
this paper in a different way. In the first step, i.e., the measurement of the ti invariant 
mass. Sec. [3], combinatorics play no role: the invariant mass can be calculated by summing 
all the final state momenta irrespective of assigning jets to top or anti-top quarks, 
or W~ bosons. In the second step, i.e. the measurement of the spin of an intermediate 
resonance, SecHJ there is in principle a 12-fold ambiguity (assuming 6-tagging) in assigning 
the correct (6-)jets to the top or anti-top quark. These ambiguities could be solved (in 
case of a very good jet energy resolution) or anyway allieviated by using the constraints 
from the top quark and W masses. In the third step, i.e., the measurement of the spin 
correlations of the top anti-top pair. Sec. O not all ambiguities can be solved: in any case 
is not possibly to uniquely identify on event-by-event basis which of the two jets come from 
the down-type quarks in the W boson decays. Experimentally, the fuUy-hadronic decay 
is difficult to trigger and extract from multi-jet backgrounds, which makes this channel 
challgenging for BSM physics studies. [78] 

The single-lepton decay channel is much more promising. The single lepton in the fi- 
nal state greatly improves the possibility for triggering on these events and extracting it 
from backgrounds compared to the fuUy-hadronic decay mode. The presence of a missing 
neutrino in the final state entails that only 5 x 3 = 15 independent measurements can be 
obtained, one short of 16 necessary. The missing information can be recovered by including 
a constraint coming from, e.g., the W boson mass (up to a two-fold quadratic ambiguity, 
which can be solved in various ways, e.g., see Ref. [31]). Using also the constraint from 
the top mass removes the ambiguity for assigning the correct 6-jet to the top quark needed 
for second and third step in our analysis. For the third step, however, it is still non-trivial 
to solve the ambiguity coming from assigning the correct jet to the down-type quark in 
the non-leptonic W decay. Several methods have been proposed, including, for example 
choosing the least energetic non-6-jet [TH [72]. Given its rate and the various reconstruc- 
tion studies and possibilities, the single-lepton channel is the most straightforward search 
channel for BSM physics in ti events. 

In the douhle-lepton decay mode there are two missing neutrino's in the final state. 
This makes the reconstruction of the full event kinematics challenging but certainly not 
impossible [73l HQ], [TH |21], [75l [761 [IZ] • There are four visible particles in the final state, 
two 5-jets and two opposite sign leptons, leading to 4 x 3 = 12 independent measurements. 
The additional four constraints from the W boson and top quark masses are just enough 
to set up a system of non-linear equations to solve for the necessary 16 variables. It can 
be shown that this system has up to eight solutions [71], of which, in general, only a few 
are physical and can be discarded or included based on their likelihood. It has to be noted 
that each solution has no further ambiguities and the event is completely reconstructed. 
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For this reason, despite the small branching ratio, this channel competes in reach with the 
single- lepton in the studies of the spin correlation studies in ti [iQl [79] . 
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